Entangling neutrons via successive scattering from a substrate 
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This letter details a simple scheme to entangle two neutrons by successive scattering from a 
macroscopic sample. In zero magnetic field the entanglement falls as the sample size increases. 
However, by applying a field and tuning the momentum of the neutrons, one can achieve a substantial 
degree of entanglement irrespective of the size of the sample. 
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Neutrons are ideal model fermions for the study of 
quantum effects. They are easily initialized to a known 
state, have well-defined degrees of freedom and a small 
interaction cross-section, which limits coupling to the en- 
vironment. Neutrons can also exhibit EPR-type correla- 
tions between distinct degrees of freedom, but so far these 
have only been observed for single particles, specifically 
between the spatial and spin parts of the neutron's wave 
function, and the neutron's energy [l[. In this letter, 
we put forward a scheme to create entangled states of 
two uncorrelated neutrons, hence opening up to further 
experimental study the fermionic analog of non-classical 
photon optics. 

The aim of our work recalls recent studies on the use 
of a macroscopic sample as both a mediator and an en- 
tanglement reservoir, as well as a vast literature on scat- 
tering entanglement mediators and sequential generation 
of entangled states (@, S 0, S, H 0, II 

and references 

therein). However, as opposed to previous proposals, our 
scheme does not require single-state preparation, manip- 
ulation, ancillary measurements or pre-existing entangle- 
ment anywhere in the system: the entanglement between 
the neutrons arises purely as a consequence of their suc- 
cessive interaction with the sample. 

We model our sample as a regular lattice of localized 
electrons, such as one might find in a ferromagnetic insu- 
lator. The sample is at zero temperature, and subject to 
a static magnetic field in the z direction. We neglect the 
nuclei, and assume the electrons interact via a transla- 
tionally invariant potential which conserves the total spin 
quantum number S z , such as a ferromagnetic Heisenberg 
exchange interaction. Disregarding boundary effects, the 
free hamiltonian of the sample then reads: 



(1) 



with J > and B z > 0, where (ij) indicates the sum 
over nearest-neighbouring pairs, J is the magnitude of 
the exchange coupling constant, B z is the strength of 
the field, TV is the number of spins in the sample, er l = 
(cr* , tjy, crj) , and the o l a are the Pauli spin matrices for 
spin i. We use the ket |0) (|f)) to denote the spin-up 
(down) eigenstate of a z , and indicate with |j) a state in 



which all the spins in the sample except that at site j are 
in state |0). 

The protocol begins by initializing the sample to a 
known pure state, chosen such that the overall state of 
the neutrons and the sample contains at most two spin 
flips. We work with sample state \ipf) = -^YlJ=i |j)> 
which is a uniform superposition of |j)-states and thus 
contains a single spin flip. Similar results can be ob- 
tained by choosing = |000..0) provided one adjusts 
the neutron polarization accordingly. Hence, we do not 
rely on \ipi) being entangled, which distinguishes our pro- 
tocol from an entanglement extraction scheme (cfr. [2]). 
Both initial states can be prepared by applying an elec- 
tron paramagnetic resonance pulse to the sample once it 
has relaxed to its ground state. In the absence of dissi- 
pative processes, the states \ipf' B ) are stable, as they are 
eigenstates of Hq. In practice, their lifetimes will be lim- 
ited by the classical spin relaxation time T\ of the system. 
However, the success of the protocol will not be affected 
provided the entangling part of the scheme is completed 
within time T\. 

After a certain period Tf <C T\ of free evolution under 
the effect of (Jf}, the sample is irradiated with a beam 
of ultra-cold neutrons (UCNs) with momentum k z z, pre- 
pared in |0) for state \^f), or a|0) +/3|1) for state j^f ), 
with a = f3 = in zero field, and a = Q, /3 = 1 above 
a threshold field B t , defined below. We assume the in- 
tensity of the source is low so only one neutron at a time 
scatters from the sample, with an arbitrary time delay 
between scatterings. This is reasonable because neutrons 
are weakly interacting particles, so for a typical flux the 
probability of many neutrons scattering at once is small. 
We model the scattering events as finite-time interactions 
of the neutrons with a composite hamiltonian of the form: 



iQR, 



(2) 



where the first term is a spin-independent potential ex- 
tending over the finite volume D 3 of the sample, and the 
second, spin-dependent term arises from the magnetic 
dipole interaction between the neutrons and the sample 
fl3j |. Here, A is a coupling strength, er™ is the spin of 
neutron m, Q is the direction of the neutron scattering 
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wavevector Q = ki — kf , and cr; is the spin of the elec- 
tron at position R; in the sample. The value of A is 
determined by A = ~g n fj, n g e /j,B^oD~ 3 , where g n is the 
neutron g-factor, fj, n is the nuclear magneton, g e is the 
electron g-factor, D 3 = Nclq and ao is the lattice con- 
stant. For the results shown, we define our energy and 
time units such that h = fig = fie = m e = 1, and set 
B t = 0.1XN. 

We assume the sample 'sees' the neutron as a de- 
localized wavepacket rather than a localized particle. As 
a result, the neutron can be thought to interact with all 
the spins in the sample simultaneously, which further dis- 
tinguishes our scheme from an entanglement extraction 
protocol (FIG.QJ. For this assumption to hold, the inci- 
dent and scattered neutron waves must be in phase over 
the whole sample, hence the neutron coherence volume 
must be of order D 3 . To maximize the scattered neutron 
flux we choose Q to be a reciprocal lattice vector, so that 
the phase factor in equation reduces to unity. 



where a am is the a-component of the spin of the inter- 
acting neutron, which we label with m, and the identity 
operation on the non-interacting neutron is understood. 
This hamiltonian gives rise to an exchange-type coupling 
between the neutron and the sample. 

Let us suppose the first neutron arrives at the sam- 
ple at time t = Tf and remains coupled to it for a finite 
time t. At time t = Tf + r it departs, and the sample 
undergoes a second period of free evolution t'j , which is 
followed by the scattering of the second neutron. We as- 
sume the durations of the scattering events are equal, so 
that the second neutron also interacts with the sample for 
a time r. This is realistic, as r is determined by the scat- 
tering process. The sequential nature of the interactions 
is assured by the fact that [^1,^2] ^ 0. This condition 
must hold if any entanglement is to be produced (l4| . 




B, 

FIG. 1: A one-dimensional snapshot of the scattering geometry 
at some fixed moment t > Tf. The sample sees the neutron as a 
wavepacket, the peak of which travels with momentum k z z. 

Let us detect neutrons scattered in the forward direc- 
tion, for which Q = 0. To account for the rapid variation 
of Q in the vicinity of Q = 0, we replace the double 
cross product in equation (|2J) with its average value over 
a small sphere of radius e <C R^ 1 , centred on Q = 0. 
One finds this average to be proportional to er;. If we 
absorb this proportionality constant into the value of A, 
the hamiltonian Ti. m can be approximated by the follow- 
ing expression: 

/ N N N \ 

T~Lm = Vo + A a xm ■ 2^, a x + a v m ' a y CTzm ' °* ) ' 
\ j=i j=i i=i / 

(3) 



The initial state of the system as a whole can be writ- 
ten as \4>o) = |0)2|0)i|V'i), where the subscripts refer 
to the neutron indices and is defined above. As 
\ijjo) is pure, the scattered state \ipf) can be obtained 
by straightforward time-evolution using the canonical 
operator U(H,t) = exp [— iHt], such that \ipf) = 
U(H2,T)U(Hi,T)\tpo). This is a slightly atypical way 
of treating a scattering problem, usually solved within 
the 5-matrix formalism which removes any explicit time 
dependence [15(. However, one can show that the two 
methods agree to first order in Vq and A, provided r = 
D\k z \~ l . This condition supplies us with a physical in- 
terpretation of the quantum-mechanical time parameter: 
in the present context, it is the classical time a neutron of 
momentum k z z would take to travel a distance D. This 
can be tuned simply by adjusting the neutron momen- 
tum. We commit a full analysis of this result to a future 
publication. 

We quantify the entanglement between the neutrons 
using the concurrence as defined by Wooters et al. 
and assume throughout this paper that J — j. The 
choice of J is entirely arbitrary, as it has no bearing on 
the behaviour of the concurrence. For illustrative pur- 
poses we set A = 1. Changing this value would simply 
re-scale our energy and time units. 



Through some algebra, it can be shown that the concurrence between the neutrons for initial state \4>f} is given by 
the following expression (we found no closed form for \ipf)): 

8V2NX 2 sin 2 4>t 



C(X,N,B.,t)= ~ v Wi =- V [0 2 + <P(X - AiV + B z ) - [J>* - 4A 2 7V sin 2 0r] , (4) 

(r[B z + A(l — Nj + <p\ 



r 



with: In zero field, the concurrence shows regular oscillations as 

^ = 4>{\,N,B Z ) = VB z 2 -2B,A(iV-l) + A 2 (iV + l) 2 . (5) 
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a function of r (FIG. [2]). The period of these oscillations 
is determined by the energy splitting of the eigenstates 
of Ti. m which correspond to the spin-flip being shared 
between the interacting neutron and the sample. For 
initial state \tpf) these oscillations persist at finite fields, 
but for \ipf) the behaviour of the concurrence is more 
complex. 




FIG. 2: The evolution of the concurrence as a function of the 
interaction time r for N = 4 and B z = 0. The green dashed 
curve relates to initial state \il)t) while the blue solid curve 
relates to initial state \ipf}- 

At B z = 0, the peak value of the concurrence C p falls 
roughly as TV -1 . However, this effect can be countered 
by switching on the field. For all N, C p is improved by the 
field, provided B z does not exceed a rough upper limit of 
2AiV. By maximizing Equation ([4} with respect to B z , 
one can show that the optimal field strength is in fact 
B* = \(N - 1) (FIG. E}. Equation (g|) then reads: 



C(N, X(N - 1), r) = 2| cos (2XVNt)\ sin (2XVNt). (6) 



Hence, we calculate the optimal interaction time: 



4AViV 



(J) 



For all N, the maximum concurrence is then C p — 0.77 
(FIG. [3]) ■ The corresponding state of the neutrons has 
the form \ipf' n ) = /i|00) + u\01) + £|10), with '- |2 
0.11, \v\ 2 ~ 0.67, |f| 2 w 0.22, £ w V6exp r " 
v / 2exp -1 ?'. 



Inl- 
and - 



The behaviour of the concurrence is independent of the 
period of free evolution between scatterings. This prop- 
erty may not extend to all entanglement measures, but is 
essential for the purpose of experimental implementation 
because we have no way of tuning ri. In zero field the 
reason for the invariance is clear, as Hq can be written 
as a multiple of the identity matrix I. In non-zero field, 
the interpretation remains outstanding. Quantitatively, 
the invariance under t'j seems specific to the class of in- 
teraction Hamiltonians {H2} which conserve S z . Hence, 
we infer that U(H.2,t) acts as a 'decoding' operation, 
which extracts the signature deposited in the sample by 
the first neutron. It follows that for the neutrons to be- 
come correlated the second neutron must scatter before 
this 'signature' disappears, i.e. within the spin coherence 
time Ti of the material. 
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FIG. 3: (Colour online) The dependence of the peak concur- 
rence C p on the magnitude of the field B z for initial state 
\ipt), with A = 1. Qualitatively similar results are found for 
initial state \tpf). 



To verify the neutrons have become entangled we use 
a witness operator W [l6|, [ItJ ■ It has been shown that 
if the state has some overlap with a 'target' state 
of the form \ip) = a\01) + 0\1O), there exists an optimal 
witness W op t which can be decomposed as follows [18| |: 



Wopt = a\z + z + )(z + z + \ +(3 2 \z z )(z z + a/3 (\x + x + ) {x + x + \ + \x~x ){x x \-\y + y ){y + y \-\y y + ){y y + \) , (8) 

I 



where |x ± ), [y^) and |z ± ) are the spin-up and down 
eigenstates of the Pauli matrices a x , a y and a Zl respec- 
tively. Such a witness could be measured with as few as 
three device settings provided one could detect all outgo- 
ing neutrons and measure each component of their spin. 
This might be achieved by performing a Stern-Gerlach 
experiment in an arbitrary direction, which is conceptu- 
ally possible though challenging from a technical view- 



point. 

Finally, we examine the major sources of experimental 
uncertainty, such as errors in calibrating the magnetic 
field or the interaction time. Let us set a lower limit 
of C p = 0.7 for the peak concurrence, and assume we 
operate either at optimal field or at optimal time. The 
allowed spread in r and B z can then be approximated by 
the relations At w lO^A^A^ and AB Z « Xy/N. The 
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fractional uncertainty in r is therefore independent of N, 
whereas the fractional uncertainty in B z is roughly pro- 
portional to N~*. These relations yield stringent but 
not unsurmountable experimental requirements, given 
the precision to which neutron velocities and static mag- 
netic fields can be calibrated [2(| HH H2] ■ 

Our proposal has an optical analogue in previous work 
by Haroche et al. |23[ on entangling pairs of atoms by 
exchange of a single photon in a high-Q cavity. How- 
ever, we underline two important differences. Firstly, we 
assume both neutrons are prepared in the same state, 
contrary to the requirement in [231 ] that the first atom be 
excited and the second be in its ground state. Secondly, 
we assume the interaction time for both neutrons is the 
same. These alterations render our proposal a realistic 
solid-state analogue of [23j, as it is currently impossible 
to prepare two successive neutrons in different spin states 
and with different momenta. 

So far, we have worked in natural units. Returning 
now to SI units, we use the technical specifications of 
the PF2 source of UCNs at the Institut Laue-Langevin 
in Grenoble to gauge some of the experimental require- 
ments of our scheme fioj . First, we estimate the required 
spin-relaxation time by imposing that T± be greater than 
the time taken by the neutrons to reach the sample. For 
UCNs with velocity v — 7 ms" 1 and a flight path of 
1CP 2 — 1 m this might require T\ as 10~ 2 — Is, which 
is achievable in materials such as phosphorus-doped sili- 
con or N@C6o [HI, EH- Second, we require that the phase 
coherence time T 2 be greater than the time between scat- 
terings. For a neutron flux F — 10 8 m s and a sample 
area of order 10 -2 m 2 we find X2 w 1 (is, also attainable 
at low temperature [111 ]. 

Next, we require the neutron coherence volume to be 
comparable to the size of the sample. This condition 



yields an uncertainty relation between coherence length 
and momentum along a certain direction a, such that 
Ap a « HAL^ 1 . Assuming our sample were, say, 10 cm 
long, enforcing this condition would require the neutron 
velocity to be exact to one part in 10 6 , which is challeng- 
ing but perhaps not unrealistic given recent progress in 
neutron spin-echo spectroscopy [24j. 

Finally, we address the structural properties of the 
sample and the robustness with respect to experimen- 
tal uncertainties. Given the form of A, B* and r*, one 
finds B* w 10^ 32 (^)- 3 T, t* w 10 21 (^)- 3 A^ s and 
D = vt* . Assuming an optimal field of 10~ 2 T and a 
sample 10 cm long, these relations yield a,Q = 10~ 10 m 
and N = 10 14 , which are both attainable values. The 
allowed spread in the magnetic field and the neutron ve- 
locity is then (AB Z /B* Z ) w 10~ 9 and (Av/v) 10 -1 . On 
both counts, this level of precision is within the capabil- 
ities of current experimental apparatus 0, H3, HI] • 

In conclusion, we have presented a simple scheme to 
create measurable entanglement between uncorrelated 
neutrons. An experimental realization would certainly be 
challenging, owing to the difficulty of detecting forward- 
scattered neutrons and performing arbitrary measure- 
ments on their spin. However, given the speed of progress 
in the field, such an experiment is perhaps not far beyond 
the reach of current neutron scattering facilities. 
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